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Abstract 

We present a conundrum that results from the imprecise use of notation for partial derivatives. Taking an 
example from mechanics, we show that lack of proper care in representing partial derivatives in Lagrangian and 
Hamiltonian formulations paradoxically leads to two different values for the time derivative of the canonical 
momentum. This problem also exists in other areas of physics, such as thermodynamics. 
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We generally insist that our students speak and write with precision, but they often do not. One case in point 
is that of partial differentiation, an operation that has a ubiquitous presence in many areas of physics, especially in 
areas of thermodynamics and canonical transformations. When we represent the partial differentiation of a function 
4'{xi,X2, - ■ ■ ,Xn) of several independent variables with respect to one of those variables Xk by we often do 
not underscore the constraints in this variation; i.e., the variables that are kept constant while Xk changes. As the 
following example from classical mechanics shows, this omission can frequently lead to serious misunderstanding. 

The comparative study of the Lagrangian and Hamiltonian approaches to classical mechanics is a standard topic 
in undergraduate courses in mechanics. The crux of the difference between the two approaches is that they have 
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two different sets of independent variables. This aspect is not always sufficiently appreciated by students. Here 
we present a simple paradox that challenges the reader to engage deeply, and thus develop a clearer perspective on 
each formalism, and their differences. Paradoxes provide many benefits in instruction, namely to motivate deeper 
thinking and provide a more thorough understanding of the topic. As pointed out by Welch [1] in this journal, their 
greatest benefit is to address "specific deficiencies in understanding and reasoning." We shall see that the resolution 
of our paradox lies in the careful handling of partial differentiation. Hence, one of our aims is to urge students to 
be cognizant of the roles played by various variables in the variation of a function that depends on more than one 
physical quantity. 

In undergraduate classical mechanics courses, students are generally encouraged to express the Lagrangian, L and 
Hamiltonian, "H in various coordinate systems. This not only helps them build very important expertise in choosing 
an efficient set of generalized coordinates for a system, often using symmetries, it also helps them explore how £ and 
Ti. could look very different in one coordinate system compared to another, and yet carry exactly the same physical 
information. 

Let us consider a Lagrangian of a free system; i.e., the potential ^(g) — 0. The Lagrangian is given by 

C{qA)^T{qA)-V{q)^T{qA) (1) 
The canonical momentum is defined by 

Then, the Lagrange's equation of motion ^ f| = implies 

dT 

Let us do a similar analysis using Hamiltonian formulation. Assuming the kinetic energy to be a homogeneous 
quadratic function of velocities, and since the potential V is identically zero, the Hamiltonian isT + y = T[5J|3]. 
Hamilton's equation of motion then implies, 

dU dT 

^-~^--^q 

Thus, the conundrum is that p is given by ^ in the Lagrangian formulation, and by — ^ in the Hamiltonian 



formalism; i.e., 

(5) 



dT 
P = ^ 



Lagrangian Formulation 



dT 
dq 



Hamiltonian Formulation 



This suggests that p — 0, as one would expect in a case of a free particle. Is this necessarily true? For example, a 
particle described by cylindrical coordinates (r, 0, z), the kinetic energy is given by T = ^mf^ + ^mr'^ip^ + ^mP , 
the derivative ^ 7^ 0, and hence pr ^ 0. Thus, for this case, it becomes difficult to explain how a non-zero Pr be 
equal to ^ in Lagrangian formulation and — ^ in Hamiltonan formulation. 

Similarly, for a particle described by spherical coordinates {r,9,(f>), the kinetic energy is given by T = ^mf^ + 
^mr^6^ + ^mr'^ sir? (fp' , neither the p^, nor the pg are equal to zero. So, how do we explain that pg is equal to 
^ in the Lagrangian formulation and — ^ in the Hamiltonian approach, and yet ps ^ 0? As we will show later, 
the solution for this apparent problem lies in the fact that ^ calculated in two different formulations are different. 
It is crucial to know which variables are being kept constant in these partial differentiations. For the Lagrangian 
formulation, the derivative ^ implies that T is being differentiated, keeping coordinates qj with j ^ i constant, and 
all (7j's are held constant. For the Hamiltonian formulation, the derivative ^ implies that T is being differentiated 
keeping coordinates qj with j ^ i constant, and all Pj's are held constant. Since the constraints under which these 
two derivatives are being taken are different, the derivatives themselves may be exactly equal in magnitude and 
opposite in sign, and hence may make Eq. ([5]) consistent with a non-zero p. 

Let us check this possibility with a free particle in spherical coordinates. The Lagrangian and the Hamiltonian 
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of the system are given by £ = imf^ -I- ^mr^O^ + ^mr^ sin^9(jP and T-L — ^ + + 2mr^sin^ e respectively [?]. 
From Lagrange's equation, we get 

dT 



= nir^ sin 9 cos 9 (j)^ . (6) 

{r,(^,r,^,0} 



From Hamilton's method, we get 



dT 



P0 



cosf 



{r,0,p.,p«,p^} TOr2sin3( 



(7) 



Substituting for p^ in Eq. ([7]), we see that pg is indeed equal to mr^ sin 6* cos9 (fP' in both cases. 

Thus, in spherical coordinates, Eq. ([5]) does not lead to any inconsistency for pg. We now show that this is 
indeed true in all coordinate systems, provided that the kinetic energy is a homogeneous quadratic function of the 
generalized velocities. 



Let us start with a system with zero potential □ and with iV-degrees of freedom, described by the Lagrangian 

N 
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(8) 



where the symmetric coefficients Fij are, in general, functions of A^-independent coordinates g,;. From here on, we 
will suppress the g-dependence of Fij . The canonical momenta pi are given by 

dC{qi,qi,t) 



Pi 



— Fij qj — Fij qj , 



(9) 



where summation over repeated indices is assumed. Since coordinates qi are assumed to be independent, the inverse 
of the matrix Fij exists, and allows us to invert Eq. ^ to find velocities in terms of momenta. We denote the 
elements of this inverse matrix by .. From Eq. solving for the velocities, we get 



9. = {F-'),^ P, 



(10) 



Since the Hamiltonian must be written as a function of qi^Pi and the Hamiltonian for the system is then given by 



'^{(li,Pi,t) = Piqi - C= Piii - -^Fijcuqj 



(11) 
(12) 



where we have substituted velocities from Eq. (|10p into Eq. ([TT]). Now, using the associativity of matrix multiplica- 
tion, we get 

H{q.,p..,t) = p, p, - i (F,, {F-') J p, ((F-i)^., p,) 



= p.{F-')^^P,--pu {{F-')^,Pt) 

= \P^ PO ' 



(13) 



where we have relabeled the dummy indices in the second term of the second line. Since the matrix Fij is a function 
of Qi's, the Hamiltonian has been properly represented in terms of qi, pi and t. Now, let us compute pi for this 
general case using the two different formalisms. In the Lagrangian formalism, from Eq. ([8]), we have 



Pi 



dC dT 1 



d_ 

dqi 



(14) 



•^It is worth noting that this discrepancy in the values for p returned by the two formaUsms is present even for the case with a non-zero 
potential V{r, 6, </>); it just appears more dramatic when the potential is zero. 



On the other hand, from Eq. (1131) . the Hamiltonian formahsm gives 

dn dT 1 ( d , \ 

Now, substituting for momenta from Eq. (|9]), we get 

dT 1/9 \ 



(15) 



9j 9fe 



(16) 



In Eqs. (I14p and (jl6p . we have computed the time derivatives of momentum pi from Lagrangian and Hamihonian 
formahsms respectively. Consistency of these two expressions requires that factors within square brackets of these 
two equations be equal; i.e.. 



d_ 

dqi 



3k 



d 



iF-^' 



F, 



mk 



(17) 



Multiplying the above equation by (_F . and summing over the index j, we get 



(F-) 



-F, 



jk 



(F-) 



mk 




(18) 



Thus, consistency of the two results requires 



S3 



-F, 



jk 



dqi 



f-')Jf,, = o 



(19) 



However, the left-hand side of Eq. is a perfect differential that can be written as 



dq, 



Fjk I = —S^ 



sk : 



(20) 



which is identically zero. Thus, for a system where the kinetic energy is a homogeneous quadratic function of the 
velocities, the two seemingly different expressions for p^, one derived from the Lagrangian and the other from the 
Hamiltonian formulation, are thus proved to be identical. 
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Conclusion: In this note, we present a conundrum tliat appears in coordinate systems where the kinetic energy is 
a quadratic homogeneous function of velocities and a function of generalized coordinates. The conundrum is most 
clear in the absence of an interaction, as both the Lagrangian and the Hamiltonian are equal to the kinetic energy, 
and the expressions for pi in two formalisms are |^ and — ^ respectively. The purpose of this note is to show that 
if C and are expressed correctly in their respective sets of independent variables, and the differentiations are done 
carefully, then there is no discrepancy. This type of puzzle goes a long way toward instilling a deeper understanding 
of the partial differentiations involved in these two different formulations of classical mechanics. We would like to 
stress that this problem is not localized to the area of mechanics alone. There are many examples in other fields as 
well. For example, in thermodynamics, the molar specific heat of a gas ^ could refer to either Cp or C„ depending 
on whether pressure or the volume is kept constant during the variation of temperature. Similarly, as pointed out 
by one of the referees, for an ideal gas we have dU = ^nRdT = ^pdV. Thus, while fp^l^ = |p, f^lr ^ These 
examples drawn from different areas of physics clearly show that for partial differentiation, we need to pay great deal 
of attention to the variables that are constrained to a constant value. 
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